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Введение. Проектирование мобильных машин невозможно без учета колебательных параметров их агрегатов, 
что требует разработки специализированных динамических моделей, учитывающих вероятностный характер 
этих параметров. Основная причина возникновения колебательных воздействий обусловлена неровностью 
профиля пути мобильной машины и непостоянством физико-механических характеристик почвы. Задачи, 
учитывающие эти особенности, с достаточной точностью решаются линейно, но в многомерных динамических 
системах такой подход недопустим ввиду наличия большого количества взаимных связей. Цель работы — 
сравнительный анализ эффективности существующих методик расчета статистических характеристик 
некоррелированных внешних воздействий применительно к мобильной машине, представленной в виде 
многомерной динамической системы с воздействиями, имеющими различные корреляционные связи. 
Материалы и методы. Внешние воздействия в многомерных динамических системах рассмотрены в 
матричной форме. При вычислении статистических характеристик учтены взаимные связи в матрицах 
спектральных плотностей, определены элементы главных и побочных диагоналей, учитываются 
корреляционные связи между воздействиями. Эти особенности значительно усложняют вычисления. Поэтому, 
чтобы получить матрицы некоррелированных воздействий, матрица внешних воздействий приведена к 
диагональному виду. 

Результаты исследования. Применением различных методик расчета проведено численное сравнение 
спектральных плотностей и интенсивности колебаний мобильной машины при варьировании скоростей 
движения и характера микропрофиля агрофона. Определенные характеристики спектральных плотностей и 
колебаний мобильных машин агропромышленного комплекса позволили выработать рекомендации по 
практическому применению представленных зависимостей для проектирования данной техники. 

Обсуждение и заключения. Результаты решения матрицы спектральных плотностей внешних воздействий 
различными методами представлены на графике спектральных скоростей колебаний. Анализ графических 
зависимостей показал, что идентичные результаты, вне зависимости от метода расчета, получаются только для 
машин со слабыми функциональными связями при некоррелированном внешнем воздействии агрофона. Для 
некоторых случаев установлены резонансные скорости движения машины. Графически показано влияние 
неровностей агрофона на интенсивность колебаний узлов машины и рассеивание нагрузок на агрегаты. 
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Введение. При проектировании надежной мобильной техники необходимо учитывать не только 
статические, но и динамические нагрузки. На характеристики динамических нагрузок значительное влияние 
оказывают колебания машины и ее узлов. Учет и определение колебательных параметров требуют разработки 
специализированных динамических моделей, рассматривающих эти параметры в аспекте вероятностного 
анализа. Основная причина возникновения колебательных воздействий обусловлена неровностью профиля 
пути мобильной машины и непостоянством физико-механических характеристик агрофона. С достаточной 
точностью задачи, учитывающие эти особенности, решаются линейно, что недопустимо при решении 
многомерных динамических систем с большим количеством взаимных связей. 

Материалы и методы. На качество технологического процесса, выполняемого мобильной 
сельскохозяйственной машиной, влияют как характеристики ее рабочих органов [1-5], так и характеристики 
ходовой системы [6,7] и микрорельеф агрофона [8, 9]. В общем случае внешние воздействия на машину 
рассматриваются как система с и функциями. Процессы изменения микронеровностей пути (или изменения 
плотности почвы) обозначим как функции 41 (Е) ...9 (Е) [9, 10]. 


Если статистическое описание внешних воздействий допускает стохастическое ортогональное 
представление, то справедливо уравнение: 


ч() =т [4.0] + [0 (9)-$. (вв) 45, 


где О;(@) — система случайных функцийс; ф; (1,0) — система неслучайных функций { и параметра 02; т [4;(1)] 
— математическое ожидание. 
Мобильные машины АПК — это системы с запаздывающими воздействиями [2, 6, 9], которые 


характеризуются многими параметрами неровности микрорельефа поля или твердости почвы [2, 8, 9]. Эти 
параметры являются однотипными и могут быть сдвинуты во времени 10. Помимо этого, конструкция машины 
(колеса переднего и заднего моста) влияет на взаимосвязь между воздействиями. С учетом этих особенностей и 
предположением, что математическое ожидание внешних воздействий равно нулю, матрица спектральных 
плотностей внешних воздействий описывается следующим образом [11, 12]: 
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5, (©) 5, (©) р Ку (1®) 
’ (<, 15) = 5, (®)е”“* 5, (©) Сл, (16) 2 (1) 
Саб уя (15) Со (1®) и (<) 


где с. =65с», =е ^%; с =1;с, =е/^* — для переднего ведущего моста; с; =е 


—Лоц . 
5 


= =е”“%. 6, =| —ДИЯ 
заднего ведущего моста. 

Для упрощения вычисления статистических характеристик приведем матрицу к диагональному виду, 
чтобы получить матрицу некоррелированных внешних воздействий. Введем характеристическую матрицу и 
приведем матрицу (1) к виду: 

5, (<) г. ^ (75) 5, (®)е ^^ Сб, (1®) 
+ + Ро, . . 
5, (в, /0)-^(15)-Е=| 5, (®)е^*  5,(®)-^(1ю) сб, (1) |, (2) 
сб (76) — с,5„(1) 5, (6)-№(до) 
где 4 — характеристическая функция; Е — единичная матрица. 
Характеристическое уравнение представим в виде: 


а? (16) +Ь^. (1%) +с^.(1ю)+4=0. (3) 
Характеристические функции ^. ( 1®),^., (15), ^. ( 1®) определялись с использованием формул Кардано 

[13]. 
Матрица перехода к новому базису, построенная по координатам векторов ^. ( 1®) (при 2 =1, 2, 3), 


определенных из уравнения (3), имеет вид: 








Ти 1 ТП 
Т(/2)= "> п» Пы|- (4) 
Тз 13 Пз 
В этом случае обратная матрица имеет вид: 
ч( и) 70°) 
р (15) —- “(6 > (5) 
где Г(/%) — присоединенная матрица; |Т ( 9) — определитель матрицы. 
Искомая матрица определялась по выражению: 
5, (18) =Т ' (1®)-5, (в, /2)-Т (15) = 41а |5, (16)-5> (16) 5, (15) , (6) 








где 51, 52, 93 — спектральные плотности некоррелированных воздействий функций комплексного переменного. 
Нахождение диагональной матрицы (6) упрощается решением плоской задачи. В этом случае внешние 
воздействия представим в виде: 














| 5, (©) 5, (®)е 
5, (в, /ю)=| и, } 7 
9 (о ©) К (=) ее“ в. (©) ( ) 
Тогда характеристическая матрица определится выражением: 
5, (®)-^(°) 5, (®)е"* 
5 —^(<).Е=|" ь | 8 
Ч (с, 5) (<) о (©) е®% 8. (©) —^.(®) ( ) 














В свою очередь, характеристическое уравнение представим как: 

^^ (©)-25, (®)-^=0. (9) 
Откуда корни уравнения равны: /, (©) = 0,^,, (©) =25, (©) : 
При корне 4, (®)=0, система уравнений для определения координат собственных векторов имеет 


следующий вид: 
5, (®) пи, (12) +5, (®)е ^^ "Пр? (1°)=0. (10) 


ЛЮ 


Задавшись 1 ‚ находим 1|> (16) =1. 


— е 
При корне 4. (©) = 25, (©) ‚ получим систему уравнений: 
| (©) 1, (762) +5, (в) е`^* п» (25) = 
5 (в) пы (1) +5, (®)-т (76) =0 


Тогда матрица перехода к новому базису принимает следующий Вид: 


(О 
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— 64 Го 
е 71° е/ о 


1 (22) = (12) 


—1 е?/®% 


Преобразованием выражения (8) найдем спектральные плотности некоррелированных воздействий: 
. | — 16 — 
5 (1) = 5, (®)-(+е ^^) =0. 
ил) — Л) — 
5, (15) = 5, (©). 1+е^*)=0. (13) 
Если допустить, что профиль агрофона в продольном и поперечном сечении определен 
некоррелированными случайными процессами [8,9], то можно представить третий элемент диагональной 


матрицы, как 5, (©) . В этом случае матрица (13) примет вид: 
5, (16) = ав |5, (© )-(1+е ^^); 5, (©). @+е^*); 5, (<}] (14) 


Решением матрицы (1), совместно с матрицей (14) и с учетом формул, описывающих некорреляционные 
воздействия, определялись спектральные плотности колебаний. В итоге получилось выражение: 


5. (©) = 5, (в) [Ф!, (©)+Ф› ()|=0. (15) 

Результаты исследования. На рис. | представлены графики спектральных плотностей вертикальных 

скоростей передней части молотилки комбайна семейства «Вектор». Очевидно, что динамические параметры 

машины влияют на спектр колебаний (рис. 1). Узкополосный спектр характерен для условий близкого 

совпадения частот собственных колебаний мостов комбайна с жаткой (рис. 1). Также на ширину спектра влияет 
запаздывание воздействий от колес мостов комбайна. 
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1 1,2 1,4 1,6 1,8 2 2,2 2,4 2,6 2,8 3 
©, -10 (с!) 
Рис. 1. Спектральные скорости колебаний при решении с использованием: 
1 — матрица (1); 2 — матрица (6); 3 — матрица (15) 


Результаты расчета интенсивности колебаний с использованием матриц (1), (6), (15) графически 
представлены на рис. 2, что подтверждает зависимость функциональных связей в машине от характера внешних 
воздействий. 
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Рис. 2. Графики интенсивности колебаний при движении машины 
по микропрофилю с различным характером неровностей: 
1 и 2 — вдоль борозды; 3 — поперек борозды поля 
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Решение систем, представляющих собой некоррелированные внешние воздействия на движущийся 
машины со слабыми функциональными связями, дают близкие результаты вне зависимости от используемых 
методов. 

Одинаковые результаты при использовании различных методов расчетов получаются только для машин 
со слабыми функциональными связями, при движении по микропрофилю поля, создающего некоррелированное 
внешнее воздействие. 

Оценка динамических характеристик комбайнов производилось по следующей формуле: 


И, = 15, (6-4, (16) 
По 


где 1, 2. — границы частот спектра. 
Интенсивность колебаний для различных условий, определенная с использованием формулы (16), 
представлена на рис. 3. 
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Рис. 3. Результаты инженерной оценки динамики комбайна 


На рис. 3 представлены зависимости перегрузок — от скорости движения комбайна. Очевидно, что 


различный спектральный состав неровностей микропрофиля агрофона влияет на интенсивность колебаний 
узлов машины. Следует отметить наличие резонансных скоростей движения в интервале от 10 до 15 км/ч. 
Полученные графические зависимости (рис. 2, 3) позволяют оценить разброс нагрузок на агрегаты. Затененная 
область иллюстрирует их рассеивание и позволяет оценить коэффициенты вариации (2). В таблице приводятся 
параметры распределения нагрузок при средних скоростях движения машин [2]. 
Таблица 1 
Динамические параметры нагрузок для мобильных машин АПК 























Тип Скорость 2. 

Режим работы ы 
машины движения м/с 8 

й 4 
Комбайн «Вектор» И - — __ 
транспортный 2,781 0,835 
МТА трактор К-700А, тележка рабочий 3,352 0,474 
2 ПТС-4 транспортный 8,323 0,510 
МТА трактор рабочий 1,247 0,364 
Т-150А, комбайн ККУ транспортный 3,527 0,382 
МТА Т-150А 

и. 29 Рабочий (пашня) 0,752 0,346 




















Обсуждение результатов. Анализ результатов расчетов спектральных плотностей колебаний по 
формулам (1), (6), (15) графически представлен на рис. 2, 3. Из графических зависимостей видно, что спектр 
воздействий нагрузок и их динамика оказывают влияние на погрешность вычислений. 
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Расчет динамических многомерных систем и оценка колебаний машин автоматизирован и производится 
с использованием специализированного программного продукта’. Программа применяется для выбора 
колебательных параметров контуров и агрегатов мобильных сельскохозяйственных машин; расчета плавности 
хода зерноуборочных комбайнов; для проектирования деталей узлов машин [14]; оценки условий труда 
механизатора. Перед использованием программного продукта необходимо предварительно выбрать 
динамическую модель; задаться конструктивными параметрами машины и ее агрегатов; определить массово- 
геометрические и упруго-диссипативные характеристики системы. 

Выводы. Выражения (1), (6), (15) приемлемы для проектировочных расчетов спектральных плотностей 
воздействий на узлы машины. Использование формулы (6) недопустимо при проектировании машин с режимом 
движения по неровностям агрофона с низкочастотным спектральным составом и для машин с низкой 
корреляцией колебаний ходовой и управляемой части. 
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